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Abstract The search for a point set configurations of the R 3 space which contains the smallest value of the 
Euclidean Steiner Ratio is almost finished. In the present work we introduce some analytical meth- 
ods which aim to support a famous conjecture of Discrete Mathematics literature. The relations of 
' this problem with that of molecular architecture in terms of motivation as well as application of the 

\ formulae obtained is also emphasized. 
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^ ! 1. Introduction 

One of motivations for modelling the macromolecular structure is the careful observation of 
the protein data banks. In the present stage of macromolecular evolution, the site of the atoms 
in a biomacromolecule can be modelled by Steiner points of an usual Steiner Problem. This 
is the Steiner Tree modelling with three edges converging at a Steiner Point and each pair 
of edges making angles of 120° therein. This looks as a last stage in the Nature's plans for 
^ ■ modelling these structures since it can be shown that this tree conformation is the most stable 

one [1]. After an exhaustive series of computational experiments by working with evenly 
spaced given points along a right circular helix of unit radius or, 

Pj = (ca&(ju), s'm(juj), aju); < j < n - 1 (1) 

where u> is the angular distance between consecutive points and 2ira is the helix pitch, we 
get a sequence of Steiner points in a helix of smallest radius and the same pitch 2na, i. e., 
belonging to the same helicoidal surface. All the experimental results led to write the formula 
for angular distance of consecutive Steiner points along the internal helix as [2] 



Uk = arctan I — I + 2tt 

Xk 



Zk 

2ira 



+ 7T 



m 
~2 



1 < k < n - 2 (2) 



and ujk+i — ujk ~ buj + c, with b « 1, c « 0. The squared brackets stand for the greatest integer 
value and m = 1, 2, 3, 4, according to the quadrant of the angle arctan 



Uk. 

We have worked with a modified version of a famous algorithm [3]. This reduced search 
space version of the algorithm was developed in 1996 and is available at http://www.biomat.org 
/apollonius. 

From eq. (2), we write the Ansatz for the Steiner Points Sk(xk, yk,Zk) as 

Sk(r(uj,a)cos(ku}),r(uj,a)sin(ku}),akui), 1 < k < n — 2. (3) 
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By forming a Steiner tree with points Pj and St with a 3-Sausage's topology [4], we can 
obtain from the condition of edges intersecting at 120°, an expression for the radius function 

r(uj, a) or: 

r(m,a)= - = (4) 

where 

Ai = l- 2cosw. (5) 

The restriction to full Steiner trees can be obtained from eq. (1) alone. We then have for the 
angle of consecutive edges formed by the given points 

4<co S ^, a ) = -l + 2(a ^ 2 + + 4 ^ 4i) . (6) 

2. The generalization of the formulae to sequences of 
non-consecutive points 

The generalization of the formulae derived in section 1 to non-consecutive points [5] along 
a right circular helix is done by taking into consideration the possibility of skipping points 
systematically to define subsequences. Let us consider the subsequences of fixed and Steiner 
points respectively. 

( P i) m ,tf max '■ P 3i P J+m, P j+2m, ■■; P j+ f max . m ( 7 ) 
(£fc)m, : Sk, Sk+m, S k+2m , •Sfc+Z^.m (8) 



with < j < m — 1 < n — 1, 0<k<m-l<n-2 and 



I 3 = 



n-j-1 



m 



■ r 

' "max 



n-k-2 



m 



(9) 



where (m — 1) is the number of skipped points and the square brackets stand for the greatest 
integer value. It is worth to say that the sequences (1) and (3) correspond to (Po)i,n-i an d 
(£1)1,71-2 respectively. 

The n points of the helical point set are then grouped into m subsequences and we consider 
a new sequence of n given points which is written as 

TO— 1 

^ U^„,/;..,. do) 
3=0 

Analogously, a new sequence for the union of the subsequences of the Steiner points is 
introduced in the form 

TO— 1 

S *= U Wm,!^ (ID 
k=l 

If the points Pj+i m , Sk+im are evenly spaced along the helices, their coordinates are given 
analogously to eqs. (1) and (3) of section 1. We can write, 

Pj+im (cos(j + lm)oj, sin(j + lm)uj, a(j + lm)uj) (12) 

Sk+im {r m (uj, a) cos(/c + lm)uj, r m (uj, a) s'm(k + lm)uj, a(k + lm)u) . (13) 
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We can now form Steiner trees with the points Pj+i m and Sk+i m with the 3-Sausage's topol- 
ogy. We get from the condition of an angle of 120° between each pair of edges meeting at a 
Steiner point, the expression of radius r m (uj, a), 

r m (w, a = 14) 

where 

A m = 1 — 2 cos(mu) . (15) 
The restriction to full Steiner trees is now obtained from the points Pj+i m only and we have 

-\ < ooflM^a) = "I + 2(m 2^+ m i ) +A m ) (16) 



3. A Proposal for a Steiner Ratio Function 

After a straightforward but tedious calculation, we get that the Euclidean lengths of the m- 
spanning tree and of the m-Steiner tree for a great number of fixed points, n » 1, are given 

by 

Vsp(u,a) = n\Jm 2 o?J 2 + 1 + A m (17) 



1st = n ^1 + maco^j TT ^-j . (18) 
The usual prescription for the Steiner Ratio Function (SRF) lead us to write 



min(l + mawJ — — ) 

p{u, a) = m V 1 + Am (19) 

miny m 2 a 2 u> 2 + 1 + A m 

m 

The "min" in eq. (19) should be understood in the sense of a piecewise function formed by 

m 

the functions corresponding to the values m = 1, 2, n — 1. 

The restriction to full trees is applied by remembering that for 3-dimensional macromolec- 
ular structure there is not tree built from partial full trees. The tree which represents the 
"scaffold" of the structure is itself a full tree or it is completely degenerate. It is seen that the 
surfaces corresponding to eq. (16) for m > 2 violate the inequality there for a large part of 
their domain. The surface for m = 1 has the largest feasible domain. Our proposal for the SRF 
function of a helical point set should be written then as 

P(",a) = . V + 1 =■ (20) 

min v m 2 ot 2 oj 2 + 1 + A rn 

m 

We can also suppose necessary bounds on p(u, a), or 



^<p(^,a)<l (21) 



where the first inequality stands for the Graham-Hwangs's greatest lower bound [6] for the 
Euclidean Steiner Ratio. 
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The corresponding ^-region is given by 

arccos(l/4) < ui < 2vr - arccos(l/4). (22) 

For the region defined by eqs. (21) and (22), we can define the unconstrained optimization 
problem of eq. (20) in the form 

p(uj, a) = max p m (uj, a) (23) 

m 

where the surfaces p m (u, a) are 



Pm{UJ ' a) = VmW + l + A m (24) 

The first three surfaces meet at a possible global minimum of the surface (20) or (23). Fur- 
thermore, in this point, beyond the 3-Sausage's topology for the Steiner tree structure we also 
have a 3-Sausage's configuration for the fixed points or the familiar set of vertices of regular 
tetrahedra glued together at common faces [7]. This non trivial solution can be written 

oj r = vr - arccos(2/3); a R = V30 [9(vr - arccos(2/3))] (25) 

and 

p{uj R , a R ) = ^(3^3 + y/7) = 0.78419037337... (26) 

which is the value assumed by the main conjecture of authors of ref. [4] for the Steiner Ratio 
in M 3 with Euclidean distance. 

It is worth to notice that there is always a pair of values which give the same value of p, 

namely (to, a) and (2ttN - uj, ^w^j) > N £ Z - This means that 



uj = n + arccos(2/3), a = v30 [9(ir + arccos(2/3))] 
correspond to the same p-value, eq. (26). 

4. The Existence of a Global Minimum. The Weierstrass 
Theorem 

The use of the Weierstrass Theorem has as a requirement the definition of a compact domain. 
We shall define it in the diagram below where we have depicted the curves p m = 1, m = 
1,2,3,4,5. 

The compact region can be defined as 

(Hpi - [jHp k )f]{(uj,a)\arccos(l/4) < 2vr - arccos(l/4)} (27) 

k>2 

where Hp\, Hp k are the hypographs of the functions p\ = 1, p k = 1, respectively. Since the 
functions pk{w, a) are continuously decreasing from the boundaries of this compact domain 
towards its interior, the Weierstrass Theorem guarantees the existence of a global minimum 
inside it. 
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Figure 1. The compact region to be used in the Weierstrass Theorem u)i — arccos(l/4), ujf — 2n — arccos(l/4) 
according the Graham-Hwang's lower bound. 

5. Concluding Remarks 

This work is one of the last steps in our effort to proof the Smith+Smith Conjecture. There 
are also nice results to be reported about the study of configurations in the neighbourhood 
of this value of the Euclidean Steiner Ratio. They are associated to a continuous description 
of molecular chirality [8, 9] and their influence in the strengthening of molecular interaction. 
Studies of this subject are now in progress and will be published elsewhere. 
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